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This  interim  scientific  report  summarizes  the  research  conducted 
under  Grant  AF-AFOSR  81-0186  during  the  period  1  August  1982  to  30  July 
1983.  The  number  in  the  bracket  [  ]  refer  to  the  items  in  the  attached 
list  of  publications. 

It  is  shown  that  using  decimation  by  D  with  an  autoregressive  model 
of  order  M  for  spectrum  estimation  yields  the  same  resolution  as  a  model 
order  MD  used  with  the  undecimated  signal,  and  that  decimation  reduces 
the  computation.  An  expression  of  the  autoregressive  spectrum  for  K 
complex  sinusoids  in  additive  white  noise  is  derived.  This  expression  is 
used  to  study  the  resolution  question  in  terms  of  the  model  order  and  the 
required  signal  to  noise  ratio  [1,2].  Extension  to  higher  dimension  is  car¬ 
ried  out  [3J,  as  well  as  extentions  to  the  maximum  likelihood  method  and 
the  use  of  autoregressive  moving  average  model  [3,4], 

The  steady  state  output  error  of  the  least  mean  square  (LMS)  adap¬ 
tive  algorithm  due  to  the  use  of  finite  precision  arithmetic  is  analyzed 
[5].  The  relationship  between  the  quantization  error  and  the  error  which 
occurs  when  adaptation  may  cease  due  to  quantization  is  also  investi¬ 
gated.  When  a  multiple  of  processors  is  used  to  implement  the  LMS  adap¬ 
tive  algorithm  in  order  to  obtain  high  throughput  rate,  the  response  is 
often  not  immediate  due  to  pipelining  and  other  factors.  A  delayed  gra¬ 
dient  estimate  LMS  algorithm  is  investigated  [6],  The  implementation  of 
adaptive  equalizers  can  be  simplified  significantly  by  employing  a  finite- 
bit  power-of-two  quantizer.  The  performance  of  such  an  adaptive  system 
is  shown  to  be  comparable  with  that  using  a  full  multiplier  [?]. 

A  common  approach  to  extrapolate  a  bandlimited  signal  in  discrete¬ 
time  is  to  solve  an  under-determined  system  of  linear  equations.  Singu¬ 
lar  value  decomposition  (SVD)  provides  a  means  for  determining  the  solu¬ 
tion  using  the  Moore-Penrose  inverse.  An  expression  for  the  mean  square 
error  is  derived  [8,9].  The  expression  indicates  that  decimation  can  be 
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applied  in  the  extrapolation  problem  to  reduce  the  high  computation  cost 
of  SVD  without  degrading  the  extrapolation.  Both  one-  and  multi¬ 
dimension  cases  are  investigated  [9], 

The  minimum  eigenvalue  and  its  associated  eigenvector  of  a  sym¬ 
metric  Toeplitz  matrix  are  of  interest  in  a  number  of  applications,  includ¬ 
ing  the  sinusoidal  in  noise  spectrum  estimation  problem.  A  fast  iterative 
algorithm  is  developed  for  this  purpose  [10].  As  in  any  iterative  scheme, 
a  good  starting  point  would  greatly  enhance  this  method.  A  one-pass 
method  for  approximating  eigenvectors  associated  with  small  eigenvalues 
is  found  [11],  which  requires  0(n2^  multiplications  where  n  is  the  size  of 
the  Toeplitz  matrix.  As  a  by  product,  we  get  reasonable  approximations 
to  all  the  eigenvalues  as  well  as  upper  and  lower  bounds  on  the  smallest 
and  largest  eigenvalues  respectively. 

A  preliminary  investigation  of  implementation  of  digital  filters 
directly  on  a  single  chip  or  a  small  number  of  chips  connected  in  an  array 
has  begun  [12].  The  emphasis  is  on  bit  serial  arithmetics  and  on  the 
operation  at  the  bit  level. 
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